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Abstract 

A finite-temperature density functional approach to describe the properties of 
parahydrogen in the Uquid-vapor coexistence region is presented. The first proposed 
functional is zero-range, where the density-gradient term is adjusted so as to reproduce 
the surface tension of the liquid- vapor interface at low temperature. The second 
functional is finite-range and, while it is fitted to reproduce bulk pH2 properties only, 
it is shown to yield surface properties in good agreement with experiments. These 
functionals are used to study the surface thickness of the liquid-vapor interface, the 
wetting transition of parahydrogen on a planar Rb model surface, and homogeneous 
cavitation in bulk liquid pH2. 
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1 Introduction 



Under normal conditions of pressure and temperature hydrogen is a gas formed by H2 
molecules. The molecule exists in two isomeric forms, para (PH2) and ortho (0H2), which 
differ in the coupling of their nuclear spins: antiparallel (J = 0) and parallel (J = 1), 
respectively. At low temperatures, near the boiling point (20.3 K), hydrogen can be catalyzed 
in nearly all parahydrogen. 

The purpose of this paper is to present a free energy density functional for liquid parahy- 
drogen. Our approach is largely inspired by density functional methods currently employed 
to study liquid helium at very low temperature .'^'^'^ Indeed, the pH2 molecule and the ^He 
atom have in common that both are spinless bosons subject to a weak van der Waals in- 
teraction. While helium remains liquid down to T = K, the stable low-temperature bulk 
phase of hydrogen, in spite of the lighter mass of the pH2 molecule, is an hpc solid. This is 
a consequence of the attractive interaction between two pH2 molecules, which is a factor of 
four stronger than the He-He one. The density functional for parahydrogen presented here 
is valid in the temperature range 14 K < T < 32 K, i.e., between the triple (13.8 K) and 
nearly the critical point (32.938 K). While the functional is still reliable in the neighborhood 
of the liquid-vapor coexistence region, the extrapolation at lower temperatures to study, for 
instance, the metastable overcooled liquid, still remains an open question. 

Ginzburg and SobyanirP have suggested that liquid parahydrogen could exhibit super- 
fluidity, analogously to liquid helium, below some critical temperature that they estimated 
to be at around 6 K. To have superfluid pH2 it is therefore necessary to bring the liquid 
below its saturated vapor pressure curve. Attempts to produce a superfluid by supercooling 
the normal liquid below the triple point have been insofar unsuccessful.'^ As a possible way 
to stabilize the hquid phase of pH2 at low temperatures several authors have considered re- 
stricted geometries in order to reduce the effective attraction between molecules. Indeed, the 
lowering of the melting point compared to the bulk liquid is a well-known and rather general 
phenomenon in clusters, see e.g. Ref. (j6j). Path integral Monte Carlo (PIMC) simulations by 
Sindzingre et a/.f^ predicted that (pH2)Ar clusters with iV = 13 and 18 molecules are super- 
fluid below about 2 K. This prediction motivated two experiment^ '^'^ and a much larger 
number of theoretical studies on small (pH2)Ar clusters, see e.g. Refs. (fTTl [T^ [T3|) and ref- 
erences therein. Concerning two-dimensional pH2 systems, PIMC calculations of films have 
found that introducing some alkali metal atoms stabilizes a liquid hydrogen phase which 
undergoes a superfluid transition below ^ 1.2 K.fHl Moreover, experiments on adsorption of 
H2 films on alkali metals substrate^^ have reported the existence of a wetting transition, 
analogous to what is known for helium. 

A density functional (DF) for liquid parahydrogen at finite temperatures could help to 
describe these phenomena. Indeed, DF methods have become increasingly popular in re- 
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cent years as useful computational tools to study the properties of classical and quantum 
inhomogeneous fluids,'^ especially for large systems for which these methods provide a good 
compromise between accuracy and computational cost, yielding results in agreement with 
experiment or with more microscopic approaches. 

A first step in that direction was made by Boninsegni and Szybisz,'^ who constructed a 
density functional to study the low temperature energetics and structural behavior of H2 
films adsorbed on Cs and Li substrates. They fitted the DF parameters by using as input 
experimental informations from the solid phase, such as the ground state equilibrium density, 
the sublimation energy, and the bulk modulus. Our approach is based on a different strategy, 
as its aim is to provide a DF that can be reliably used in the hquid- vapor coexistence region. 

This paper is organized as follows. In Sec. 2 we describe the density functional for liquid 
parahydrogen in two versions, neglecting or not finite-range effects. The validity of the DF is 
assessed in Sec. 3 by comparing the calculated surface tension with experiment. The DF is 
thus used to study the wetting of a Rb planar surface by pH2, and the results are compared 
to those obtained by Quantum Monte Carlo simulations. As a final application, in Sec. 4 we 
describe homogeneous cavitation in liquid pH2. A brief summary and outlook are presented 
in Sec. 5. 

2 The density functional 

Following the procedure previously employed for liquid ^He, in particular the works of 
Guirao et a/.,'^ Dalfovo et a/.,'^and Ancilotto et a/.,'^we present first a density functional 
for bulk liquid parahydrogen. The extension to inhomogeneous liquids will be considered in 
the next section. In the most general form, the free energy of a Bose system is assumed to 
be a DF, which we write as 



where fni is the free energy density for an ideal Bose gas, and fc is the correlation energy 
density, which incorporates the dynamic correlations induced by the interaction. The number 
of pH2 molecules per unit volume is denoted by p. 

2.1 The non-interacting part 

The expression for fni{p,T) at T 7^ is given in textbooks. Following Ref. (|2T1) it can be 
written as 




(1) 



ni 



p ksT ln{z) 




95/2{z) , 



(2) 
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where 

A = ^2Txn^/mkBT (3) 
is the pH2 thermal wavelength, m is the H2 mass, and the fugacity z is defined as: 

if pA3> (73/2(1) 



if pA3< ^3/2(1) , 

where zq is the root of the equation pX^ = g^i2{z). In the above equations, gp{z) = /P. 
In the following we shall use units such that the Boltzmann constant /eb is one, so that energy 
is measured in K. In these units %^ /m = 24.0626 K A^. 



2.2 The correlation part 

Once the functional Eq. ([T]) is defined, a complete thermodynamical description of the 
bulk liquid can be achieved. In particular, the chemical potential and the pressure are given 
by 

/^(P, T) = finiip, T) + (5) 

and 

p(p,r)=p™(p,r) + p^(^^i|M^^. (6) 

In the next subsections we propose a phenomenological expression for the correlation en- 
ergy density, inspired by previous works on liquid helium, which contains a few adjustable 
parameters which are fixed so as to reproduce selected experimental properties of the liquid 
from the data collected by McCarthy et al^ or the more recent compilation by Leachman 
et al^ In fact, the latter compilation has been summarized in terms of a rather complex 
parametrization of the bulk liquid pH2 free energy density which reproduces satisfactorily 
the available experimental data between the triple and the critical points in a wide range of 
pressure values. 

We present here a simpler yet less accurate functional that has the advantage of being 
easier to generalize to finite systems (drops, liquid pH2 in confined geometries) and presents 
a more conventional separation into kinetic and potential terms. This may be of interest 
for other applications, such as computing the excitation spectrum as well as the static 
and dynamic responses of these systems.'^ Besides, the present DF could in principle be 
extrapolated to lower temperatures near T = 0, one of our medium-term goals. 
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2.2.1 Zero-range 

We consider a polynomial in powers of the density to represent the correlation part of the 
free-energy density 

Mp,T) = a,{T)p' + a2{T)p' + a,{T)p'' + a,{T)p' . (7) 

For a given value of T, there are thus four free parameters. With this form, other quantities 
such as the chemical potential p, the pressure p, and the speed of sound c, can be written as 

= fJ'ni + 2aip + 3a2p^ + ^a^p^ + 5a4p^ 
P = Pni + aip^ + 2a2p^ + Sasp^ + Aa^p^ 

mc^ = p^^ + 2aip + Qa2p^ + 12a3p^ + 20a4P^ , 
op 

where, for the sake of simplicity, we have omitted the T-dependence in the parameters. 

These parameters are fixed by the following strategy. At a given T, we determine the 
liquid- vapor equilibrium by imposing the equality between pressures and chemical potentials. 
The liquid and vapor densities are experimentally known.'^^ This provides two conditions. 
The other two conditions are the fitting of the pressure and speed of sound of the liquid phase 
at the saturated vapor pressure to the corresponding experimental values.'^ The resulting 
ai{T) parameters are given in Tabled! We shall refer to them as the DFpH2 set. 

In Fig. [1] the experimental and calculated pressure and speed of sound as a function of 
p for three values of T are compared. For each T, the first density value shown is that of 
the liquid at the liquid-vapor coexistence, where p and c coincide by construction with the 
experimental values. It is worthwhile seeing that the DFpII2 parametrization can be safely 
used for densities well above those of the liquid at the saturation vapor pressure. 

2.2.2 Finite-range 

The functional Eq. ([7]) is a zero-range one, and it cannot account for inhomogeneities or 
finite size effects induced by strong density compression. Inspired in the way this problem 
has been addressed in the case of liquid ^He,'^ we derive from Eq. ([7]) a finite-range density 
functional with the following prescriptions. On one hand, the term in p^ is replaced with 

ai(T)p2 i I dr'p(r)V(|r - r'|)p(r') , (8) 
where V represents the II2-II2 pair interaction screened at short distances 

n-) = 4e|(^)''-(^)1 if r>MT), (9) 
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and V{r) = otherwise. We have used the Lennard- Jones parameters e = 34 K and a = 3.06 
A, which are shghtly different from the standard ones (34.2 and 2.96, respectively). The value 
of h{T) is fixed by the condition 

ai(T) = i|dry(r). (10) 

The values of h{T) are given in the last column of Tabled! 

On the other hand, a coarse-grained density is introduced in the remaining terms of the 
correlation part Eq. ([7]): 

a,{T)p'+' ^ a^{T)pp' , t>2, (11) 

where 

p(r) = J dv'w{\r - r'|)p(r') (12) 

with 

and w{r) = otherwise. For the reasons explained in Sec. 3.2, we have taken h{T) = 
l.Uh{T). 



3 Surface and interface properties 

Once the DF parameters have been adjusted to reproduce the bulk behavior of liquid pH2 
at finite temperatures, the predicted surface and interface properties provide a significant 
test about its validity. We present in this section the results of a series of calculations aimed 
at showing that our proposed DF can indeed be reliably used to study finite temperature 
properties of inhomogeneous pH2 systems, such as a free-standing fluid-vapor interface or 
liquid pH2 adsorbed on a solid surface. 

We have first computed the liquid-vapor surface tension of pH2 as a function of tem- 
perature, and compared our results with the experimental data of Ref. (j^Sj). We obtain 
a liquid-vapor interface by minimizing the DF with the appropriate boundary conditions, 
corresponding to a planar liquid pH2 film in equilibrium with its own vapor. In this case the 
solution depends only on the coordinate z normal to the surface, and provides the equilib- 
rium density profile p{z). From the knowledge of p{z) one can extract values for the surface 
tension 7(T). We distinguish in the following two different approaches used to compute this 
quantity, the first one being based on the zero-range DF described in Sec. 2.2.1, and the 
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second one on the more accurate but computationally more expensive non-local, finite-range 
DF described in Sec. 2.2.2. 

3.1 The pH2 surface tension from zero-range DF calculations 

Our starting point is the zero-range free energy density functional 

f = fm + fc+/3^^^+avp?- (14) 

2m p 

This is the simplest extension of Eq. ([T]) in Sec. 2 including the first two terms in a gradient 
expansion to account for density inhomogeneities. We make here the usual choicd^ /3 = 1/4 
for the kinetic energy term, while we fix ^ to the value ^ = 19306.74 K A^, which allows to 
reproduce the experimental surface tension at T = 14 

With such a form for the free energy density, the surface tension turns out to b^^ 

7(T) = 2 r dp [/(p) - /(p.) - p(p - p.)]i/2 f- + , (15) 

Pu pv being the asymptotic (constant) density values on the liquid and vapor side of the 
interface, respectively. 

We show in Fig. |2]our results for the pH2 surface tension. The open squares correspond to 
the calculated points, while the dotted line shows the experimental result from Ref. (I22p . It 
appears that, once the DF parameter ^ is fixed to reproduce the value at a given temperature, 
say T = 14 K, the predictions for 7(T) are quite accurate in the whole range of temperatures 
between the triple and near the critical point. 

3.2 The pH2 surface tension from finite-range DF calculations 

We use here the finite-range description of the free energy density functional of pH2, as 
discussed in Section 2.2.2. Concerning the coarse-grained density defined in Eq. f[T^ . the 
empirical prescription h{T) = 1.14 h{T) is the one which gives optimum agreement between 
the calculated and measured values for 7(T). 

We use a slab geometry in our calculations, with periodic boundary conditions whereby 
a thick planar liquid film is in contact with a vapor region. Our calculations, which are 
effectively one-dimensional, are performed within a region of length Zm- Rather large values 
of Zm and of the film thickness tf (no to be confused with the thickness of the liquid- vapor 
interface) are necessary for well converged calculations (typically Zm ~ 200 A and ~ 150 
A), in order to let the system spontaneously reach (i) the bulk liquid density in the interior 
of the film and (ii) the bulk vapor density in the region outside the film. 
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The equilibrium density profile p{z) of the fluid at the interface is determined from the 
variational minimization of the grand potential, leading to the following Euler-Lagrange 
(EL) equation: 



2m dz^ 



+ u[p{z)]\^) = ^Ji^), (16) 



where the effective potential U is defined as the functional derivative of the free energy 
density^ and the value of the chemical potential /i is fixed by the normalization condition 
to the total number of H2 molecules per unit area, / dzp{z) = N/A. 

In Fig. [3] we show some selected equilibrium density profiles around the liquid-vapor 
interface calculated using the finite-range DF. As expected, the width of the interface in- 
creases monotonically with temperature. Density profiles for H2 have been obtained in Ref. 
(l26|l using the Silvera-Goldman potentiaP^ within a molecular simulation of the liquid- vapor 
interface. However, the rather thin liquid slab used there is probably insufficient to prop- 
erly describe the liquid-vapor interface. For this reason we do not attempt a quantitative 
comparison between our calculated profiles and those reported in Ref. ( |26|1 . 

The liquid- vapor surface tension 7 can be directly computed from the definition j = {Q + 
pV)/A once the equilibrium density profiles have been obtained. In the previous expression, 
Q = F — pN is the grand potential, p is the saturated vapor pressure at temperature T, V 
is the volume of the system, and A is the surface area. For the one- dimensional problem 
considered here, one has: 
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2 



f[p{z)]dz-pl p{z)dz+pzr, 
Jo 



(17) 



The pressure p can be conveniently calculated from the definition p = [pp^ — f{pv)] = 
[ppi — f{pi)]- A factor 1/2 appears in the previous equation to account for the two free 
surfaces delimiting the liquid film in our slab geometry calculations. The surface tension 
obtained from the finite-range DF is also shown in Fig. [2] as full squares. The overall 
agreement with the available experimental dataP^ is quite satisfactory. 

We show in Fig. Hlthe results obtained for the thickness of the liquid- vapor interface. We 
define the thickness as the distance across the interface between the points at which the fluid 
density changes from 0.9 to 0.1 times that of the liquid at saturation for the corresponding 
p and T values. Open squares represent the values obtained using the zero-range DF, and 
full squares those obtained using the finite-range DF. The latter are our prediction for the 
thickness of the interface. Similarly to what occurs for liquid ^He,'^'^the thickness of the 
liquid-vapor interface obtained from the finite-range DF at low temperatures is about 2 A 
smaller than that obtained from the zero-range DF. 
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3.3 Wetting properties: the pH2/Rb system 

The existence and nature of the wetting transitioii^ at a temperature above the triple 
point Tf (and below the critical temperature) is known to be the result of a delicate balance 
between the interatomic potential acting among the atoms in the liquid and the atom- 
substrate interaction potential. When the latter dominates, the fluid tends to wet the 
surface. 

When a first order transition from partial to complete wetting occurs at a temperature 
above the triple-point (and below the bulk critical temperature), then a locus of first-order 
surface phase transitions must extend away from the liquid- vapor coexistence curve, on the 
vapor side. At T < the thickness of the adsorbed liquid film increases continuously with 
pressure, but the film remains microscopically thin up to the coexistence pressure psat{T), 
and becomes infinitely (macroscopically) thick just above it. At temperatures < T < T^^, 
being the prewetting critical temperature, the thin film grows as the pressure is increased 
until a transition pressure Ptr{T) < Psat{T) is reached. At this pressure a thin film is in 
equilibrium with a thicker one, and a jump in coverage occurs as p increases through ptr{T). 
At still higher pressures this first order transition is followed by a continuous growth of the 
thick film which becomes infinitely thick at Psat{T)^^ 

Wetting transitions have been observed for quantum fluids (helium and hydrogen) on 
alkali metal surfaces. In particular, H2 on Rb surface shows a wetting transition.'^Sl This 
system has been addressed by finite temperature quantum simulations in Ref. ( 130|) . where 
a first-order wetting transition was indeed found to occur. The value found for the wetting 
temperature, Tw = 25.5±0.5 K was higher than the experimental value = 19.1 K.f^^The 
disagreement has been attributed to the fact that the interaction potential used to simulate 
the fluid-surface interaction was probably too weakly attractive. We are not concerned here 
in reproducing the experimental value for Tw but rather to compare the prediction of our 
finite-range DF approach against an accurate benchmark, i.e. the virtually exact Quantum 
Monte Carlo results reported in Ref. ( l3Up . The external potential Vs{z) {z is the coordinate 
normal to the surface plane) used in Ref. ( 1511]) to describe the H2/Rb surface interaction is a 
semi-a5 initio one.'^ In order to make a proper comparison with the results of Ref. (130]) we 
have used the same atom- surf ace potential employed there, and calculate within finite-range 
DF theory the wetting temperature for the pH2/Rb system. 

The free energy finite-range functional for pH2 described in Sec. 2.2.2 must be augmented 
by a term describing the interaction (per unit surface) of the pH2 fluid with the Rb surface, 
/ dzp{z)Vs{z), and the equilibrium density profile p{z) of the fluid adsorbed on the surface 
is determined by direct minimization of the density functional with respect to density varia- 
tions, which yields an equation similar to Eq. (1T6]) with an extra potential term arising from 
Vs{z). Details on the procedure can be found e.g., in Ref. ( 13 2 p and references therein. 
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We show in Fig. \5\ a set of calculated density profiles for different temperatures. We fix 
in our calculations the pressure (equivalently, the chemical potential) just below the bulk 
liquid saturation pressure, and increase the temperature. The abrupt (first order) change in 
the amount of liquid adsorbed on the surface shown in Fig. [5] when T is just above T = 28 
K is a clear signature of the pre-wetting transition very close to saturation. We thus take as 
our prediction for T^y the value = 28.5 ± 0.5 K. This value compares well with the one 
found by Quantum Monte Carlo calculations.'^ 

4 Cavitation 

Another application of the density functionals described in Sec. 2 is the study of thermal 
homogeneous cavitation in liquid pH2. It can be straightforwardly extended to the case of 
nucleation of pH2 droplets in a supersaturated parahydrogen vapor. From the experience 
gathered in similar ^He studies, it is enough to use the zero-range DF for this application. 
We refer the reader e.g., to Refs. [Ml ESI ESI EZ]) and references therein for a general 
discussion of similar work carried out for liquid helium. 

It is well know that phase transitions take place in the coexistence region of different phases 
of homogeneous systems, and that it does not always happen under equilibrium conditions. 
Indeed, as the coexisting phase forms, the free energy of the system is lowered, but the 
original phase can be held in a metastable state close to the equilibrium transition point. 
The metastable state is separated from the stable state by a thermodynamic barrier that 
can be overcome by statistical fluctuations. 

Cavitation consists in the formation of bubbles in a liquid held at a pressure below the 
saturation vapor pressure at the given temperature. In this case, cavitation is a thermally 
driven process. The formation of the vapor phase from the metastable liquid proceeds 
through the formation of a critical bubble, whose expansion cannot be halted as its energy 
equals the thermodynamic barrier height. A very simple estimate of the radius Rc of this 
critical bubble is provided by the thin wall model, which consists in treating the bubble 
as a spherical cavity of sharp surface and radius R flUed with vapor at the saturated vapor 
pressure, Psvp- The barrier for creating a bubble of radius R is the result of a balance between 
surface and volume energy terms 

U{R)=A7r^R'-—R^\Ap\, (18) 

where |Ap| is the difference between the applied pressure and Psvp, and can be positive and 
negative as well. This barrier has thus a maximum Umax = 167r7^/(3| App) at a critical 
radius Rc = 2-f/\Ap\. Taking 7 = 2.96 dyne/cm at T= 14 K,'221and \ Ap\ ~ 78 bar (the pHa 
spinodal pressure at T = 14 K), one gets Rc ~ 7.5 A. 
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The cavitation rate J, i.e., the number of critical bubbles formed in the liquid per unit 
time and volume is given by 

J = Jot exp{-Umax/T) , (19) 

and the prefactor Jqt depends on the dynamics of the cavitation process. A precise knowledge 
of the prefactor is not crucial, since at given T, the exponential term in Eq. (1191) dramatically 
varies with p in the pressure range of interest, thus making the determination of ph rather 
insensitive to its actual value.'^ We shall estimate Jot as an attempt frequency per unit 
volume, namely Jqt = T/{hVc), where h is the Planck constant and Vc is the volume of the 
critical bubble. 

To determine the homogeneous cavitation pressure within DF theory, one has first to 
determine the critical bubble for p and T values corresponding to points lying within the 
liquid-vapor phase equilibrium region where bubbles may appear. This is done by solving 
the EL equation arising from the functional variation of the grand potential with respect 
to the pH2 particle density. The grand potential density is ijj{p,T) = f{p,T) — fip, where 
f{p,T) is the free energy per unit volume including surface terms arising from the bubble 
surface and kinetic energy contributions already discussed in Sec. 2.2.1, and p is the chemical 
potential of the bulk liquid. The EL equation is solved imposing the physical conditions that 
p'(r = 0) = and p(r — )• oo) = pi, where pi is the density of the metastable homogeneous 
liquid. 

The nucleation barrier Umax is obtained by subtracting the grand potential of the critical 
bubble from that of the homogeneous metastable liquid at the given p and T 

Umax = Jdr[f{p,T)-Mpi,T)-p{p-pi)], (20) 

where fi{pi,T) has no density gradient terms as it corresponds to the metastable homo- 
geneous liquid. The above expression embraces the two limiting physical situations of a 
barrierless process when the system approaches the spinodal point, and of an infinite bar- 
rier process when the system approaches the saturation line. While the latter limit is well 
reproduced by the thin wall model, the former is not.'^'^ 

Figure M shows several density profiles for critical bubbles at selected T values. For each 
temperature, we have chosen the pressure corresponding to the homogeneous cavitation 
value, see below. The radius of the critical bubble for T = 15 K is similar to that obtained 
within the thin wall model at T = 14 K. 

Figure [7] shows the barrier height of the critical bubble for selected T values as a function 
of p (solid lines). The divergence of Umax as p approaches the saturation line, as well as its 
disappearance at the spinodal line are clearly visible in the figure. 
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At a given temperature, there is an appreciable probability of cavitation occuring when p 
reaches a value such that 

l = {Vx t). Jot exp{-UmajT) , (21) 

where {V x t)e is the experimental volume times the experimental time. This defines the 
homogeneous cavitation pressure ph, an intrinsic property of the liquid that is determined 
by solving the above algebraic equation. To solve Eq. f l2T|) one needs some experimental 
information. If the experiment were carried out as for liquid ^He,!^ {V xt)f, ^ 2x 10~^^ cm^ 
s. Lacking of a better choice, this is the value we shall adopt here for parahydrogen. 

Figure [8] shows the homogeneous cavitation pressure as a function of T. Also shown is the 
spinodal line Psp{T) and the saturation vapor pressure line Psvp{T)- The three lines merge 
at the critical point. It can be seen that ph is negative up to T ~ 29 K, meaning that a 
large tensile strength may be needed to produce cavitation. Notice that, as in 
homogeneous cavitation pressure is just above the spinodal pressure. 

Finally, we mention that the dashed line in Fig. [7] represents the barrier height at p^ 
for the chosen T value. It can be seen that these heights are in the 500-1000 K range for 
temperatures 14 < T < 30 K. 

5 Summary and outlook 

In this work we have proposed a density functional that can be used to study the properties 
of parahydrogen in the liquid-vapor phase-diagram region and its neighborhood. We have 
used its zero-range version to address thermal homogeneous cavitation in the metastable 
liquid, and its finite-range version to predict the thickness of the liquid-vapor interface. 
We have shown that this finite-range DF predicts the surface tension of the liquid-vapor 
interface and the wetting properties of pH2 on planar Rb surfaces in satisfactory agreement 
with experiments and microscopic calculations, respectively. 

Besides thermal cavitation, the existence of a finite temperature, finite-range DF makes it 
possible to improve the theoretical description of phenomena that hardly could be described 
using more fundamental approaches, as the study of electron bubbles in pH2.'^'^ Large 
pure or doped parahydrogen drops at finite temperature could also be addressed provided 
the impurity-pH2 interaction is known, as well as extending the study of wetting to other 
substrates and/or geometries. 

The present DF may be of interest for other applications, such as the calculations of the 
excitation spectra of pure or doped pH2 drops, as well as their static or dynamic responses. 
The extrapolation of the DF to temperatures near T = still remains an open question. We 
plan to study some of these phenomena in the future. 
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TABLE 1: DFpH2 set of parameters aj(T) (in K A^*) entering the correlation energy 
density functional fc- The parameter h{T) (in A) related to the finite-range DF is given in 
the last column. 



T 


ai/lO^ 




as/lO^ 


03/10^ 


04/10^ 


h 


14 


-5.39621105 


2 


25023742 


-1.20053135 


2.79594232 


2 


98886735 


15 


-5.28361552 


2 


15103496 


-1.14667549 


2.70797035 


2 


93284836 


16 


-5.18646074 


2 


07504193 


-1.10362936 


2.63656153 


2 


90255249 


17 


-5.10162915 


2 


01565556 


-1.06770278 


2.57547938 


2 


88158016 


18 


-5.02544156 


1 


96727851 


-1.03616377 


2.52033226 


2 


86543169 


19 


-4.95500617 


1 


92571884 


-1.00701996 


2.46790257 


2 


85212913 


20 


-4.88844024 


1 


88827147 


-0.978964063 


2.41596399 


2 


84067446 


21 


-4.82469889 


1 


85328416 


-0.951137931 


2.36290741 


2 


83053490 


22 


-4.76330170 


1 


81966306 


-0.922799978 


2.30707709 


2 


82141283 


23 


-4.70410834 


1 


78695682 


-0.893590816 


2.24744238 


2 


81313315 


24 


-4.64712668 


1 


75458070 


-0.862880207 


2.18218726 


2 


80558079 


25 


-4.59255499 


1 


72283140 


-0.830703205 


2.11075548 


2 


79868923 


26 


-4.54035421 


1 


69068884 


-0.795990139 


2.03000386 


2 


79237811 


27 


-4.49046505 


1 


65741296 


-0.757823418 


1.93673005 


2 


78657987 


28 


-4.44267216 


1 


62195349 


-0.714989528 


1.82655914 


2 


78122146 


29 


-4.39598018 


1 


58109948 


-0.664408250 


1.68988952 


2 


77615776 


30 


-4.34840634 


1 


52965533 


-0.601172851 


1.51084983 


2 


77116004 


31 


-4.29478440 


1 


45544413 


-0.514224892 


1.25423884 


2 


76570856 


32 


-4.21484337 


1 


31639928 


-0.367772398 


0.809972555 


2 


75790923 
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Figure 1. (Color online) Pressure and speed of sound: A comparison between experimental 
data (symbols) and DFpH2 results (lines) at three values of T. 
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Figure 2. Calculated surface tension. Open squares: results of the zero-range DF. Full 
squares: results of the finite-range DF. Dotted line: experimental data from Ref. (I22ll . 
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Figure 3. Density profiles around the liquid-vapor interface calculated using the finite- 
range DF. From top to bottom in the liquid regime, the profiles correspond to T = 14, 18, 
22, 26, and 30 K. The origin of the z axis is arbitrary. 
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Figure 4. Thickess of the liquid-vapor interface. Open squares: results of the zero-range 
DF. Full squares: results of the finite-range DF. The vertical hues indicate the triple and 
critical temperatures. 
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Figure 5. Density profiles of pH2 on a planar Rb surface at different temperatures. From 
bottom to top, in the asymptotic vapor regime: low density profiles, T — 22, 23, 24, 25, 26, 
27, and 28 K; high density profiles, T = 29 and 30 K. 
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Figure 6. Density profiles of critical bubbles corresponding to several T values obtained 
using the zero-range DF. For each temperature, we have chosen the pressure corresponding 
to the homogeneous cavitation value. 
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p(bar) 

Figure 7. Solid lines: barrier height Umax of critical bubbles for selected T values (K) as 
a function of p (bar) . The dashed line represents the barrier height ai p = Ph for the chosen 
value of T. 
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